SYMMETRY RESTORATION FOR ODD-MASS NUCLEI WITH A 
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In these proceedings, we report first results for particle-number and angular-momentum 
projection of self-consistently blocked triaxial one-quasiparticle HFB states for the de- 
scription of odd-j4 nuclei in the context of regularized multi-reference energy density 
functionals, using the entire model space of occupied single-particle states. The SIII pa- 
rameterization of the Skyrme energy functional and a volume-type pairing interaction 
are used. 



O !• Introduction 

^ Methods based on the self-consistent mean-field approach provide a set of fully 

microscopic theoretical tools that can be applied to all bound atomic nuclei in a 
systematic manner irrespective of their mass, N / Z ratio, and deformation.^^ Pure 
^ mean-field methods, however, have several limitations. The first one is due to the 

' determination of a wave function in an intrinsic frame of reference. Although the 

symmetry-breaking mean-field approach is a very efficient and transparent way 
to grasp the effect of correlations associated with collective modes in the limit of 
I strong correlationspP' the absence of good quantum numbers and the correspond- 

T-H ing selection rules compromises the calculation of transition moments. For example, 

broken rotational symmetry mixes states with different eigenvalues of J^, i.e. the 
. . members of a rotational band, and independent quasiparticle states of BCS type 

^ ^ used to describe pairing correlations are spread over several particle numbers. The 

second limitation concerns nuclei for which a mean-field description through a sin- 
gle configuration breaks down because several configurations with different shell 
structure are close in energy. To overcome these limitations requires what is often 
called " beyond- mean-field methods", i.e. symmetry restoration and configuration 
mixing within the Generator Coordinate Method (GCM)EEIIIII 

There are many more applications of mean-field and "beyond" methods to even- 
even nuclei than there are for odd-v4 nuclei (not to mention odd-odd ones). The 
most important reason is that the last nucleon in the odd-yl system breaks some of 
the symmetries that make such methods particularly efficient for even-even systems. 



*Based on talk presented at 18th Nuclear Physics Workshop "Maria and Pierre Curie", 2011, 
Kazimierz, Poland. 
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2. Odd-mass nuclei in self-consistent mean-field models 

Odd-mass nuclei are described as one-quasiparticle states build on a HFB vacuumP 
The most transparent way to represent such a state is in the canonical basis {a|} 

kb)=al n {U^^ + V^-^alal)\-) (1) 

where i and i are the indices of the paired conjugate states, where we use the usual 
convention i > and i < 0, and b is the index of the blocked single-particle state. 
The single-particle states and Bogoliubov matrices U and V are determined self- 
consistently. The corresponding HFB equation is obtained from variation of the 
energy with the obligatory constraints ensuring that jg^) remains a quasiparticle 
vacuum and on proton and neutron numberEEl Optional constraints are on the 
two components of the quadrupole tensor in the intrinsic major axis system of 
the triaxial shape and on the angular momentum component along one major 
axis. As done in our earlier configuration mixing calculations, we add the Lipkin- 
Nogami prescription to the HFB equations to enforce pairing correlations in all 
mean-field stateslSEi The resulting HFB equations for the blocked state are solved 
using the "two-basis method" , where in an iterative procedure the HFB Hamiltonian 
is diagonalized in the single-particle basis of eigenstates of the HF Hamiltonian in 
a coordinate-space representation.!^ 



3. Going "beyond the mean field" 

Restoring the symmetries decomposes a given SCMF state into states with good 
quantum numbers. Eigenstates of the particle-number operator N with eigenvalue 
Nq are obtained applying the particle-number projection operatoil^ 

Pn,,^^ fj d^Ne^^-^^-""-) (2) 

to the SCMF states for protons and neutrons. Eigenstates of the total angular- 
momentum operator in the laboratory frame and its z component Jz with 
eigenvalues h?J{J+ 1) and Mf, respectively, are obtained applying the operatoil^ 



2J+ 1 



167r2 7„ 
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P'iiK = =^7rV / da I dp sin(/3) / d^i Vj/j, {a, P, 7) R{a, /3, 7) , (3) 



where R = Q-^<^-h g;-tpJy ^-n^h jg ^j^g rotation operator and 'Dfjj^ a Wigner func- 
tion. Both depend on the Euler angles a, /3, and 7. -P/^^ picks the component with 
angular-momentum projection K along the intrinsic z axis from a mean-field state. 
The projected state is then obtained by summing over all K components 

+ ./ +.7 

\JMvq^,)= J2 f'iK,qh)PitK\Qb)= J2 f'iK,qb)\JMKqt) (4) 

K=-J K=-J 
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with weights f^{K,qi,) determined from a variational principle that leads to the 
so-called Hill- Wheeler-Griffin (HWG) equatioiJUH 



which for sake of simple notation we write for a Hamilton operator. In practice, 
we use an effective interaction that is provided by a multi-reference (MR) energy 
density functional (EDF) . It is the generalization of the single-reference (SR) EDF 
employed in SCMF methods to calculate the non-diagonal kernels entering the 
symmetry-restored energy and the HWG equation ([s]). It is usually postulated in 
a form that preserves analogies with the case of a Hamilton operatorEEEni How- 
ever, the MR EDF has to be regularized to remove spurious contributions to the 
energy that manifest themselves through divergences and/or finite steps when plot- 
ting the symmetry-restored energy as a function of a collective coordinate We 
have implemented such regularizatiorEMll into our codes for general configuration 
mixingf^ We use the parameterization SIlP^ of the Skyrme EDF together with a 
contact pairing interaction of volume type and a strength of 300 MeV fm^ together 
with cutoffs at 5 MeV above and below the Fermi energy!^ The Coulomb exchange 
term is neglected as done in our earlier regularized MR EDF calculationsEllHl The 
non-diagonal norm kernels entering the symmetry-restored HWG equation ^ are 
calculated directly with their sign in a technique based on a Pfaffian. 

4. An illustrative example: ^^Cr 

We have chosen ^^Cr for a first exploratory study as it exhibits several low- lying 
collective rotational bands that can be easily associated with particular blocked 
single-particle levels. A detailed study of the adjacent even-even ''^Cr that allows 
for the analysis of the change brought by the additional nucleon is also underway. 

4.1. Projected energy surfaces 

The first step of our calculation was to determine which quasiparticle gives the 
lowest energy at the Single-Reference level. It is a quasiparticle with parity minus 
and with the mean value of the projection along the deformation axis of the angular 
momentum close to | : {J^}'^ ~ 2.5^, which agrees with the experimental finding 
for the ground state that has = | . We then tried to follow this quasiparticle 
in the first sextant of the /3 — 7 plane. As we have neither good J nor good J3, 
however, it is not always clear which quasiparticle should be blocked. This ambigu- 
ity presents an additional motivation to go to the Multi-Reference level, that is to 
restore the angular momentum. As demonstrated for example by Schunck et al.^^ 
because of the time-odd components in the EDF the sextants are not equivalent, 
and a full study would require to explore at least three out of the six. Test calcu- 
lations in the other sextants indicate, however, that for this blocked quasiparticle 




K' = -.] 
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the differences after K mixing remain tiny, and do not have any importance for the 
present discussion. This might not be the case for other blocked quasiparticles, and 
certainly there will be significant differences between the sextants when cranking 
the quasiparticle states to higher spin. 



Non-projcctcd ^ 




Fig. 1. Top left: Non-projected energy surface, in the first sextant of the /3 — 7 plane, of the 
lowest one-quasiparticle state with (Js)^ ~ 2.5~ in ^^Cr. Top right: Particle-number restored 
energy surface constructed from the same one-quasiparticle states. Bottom: Particle-number and 
angular-momentum restored energy surface for the lowest J'^ = | , based again on the same 
one-quasiparticle states. Energies are normalized to the minimum of the non-projected surface. 
The deformation parameter /3 = (Q2) with R = 1.2A^^''^ fm and 7 are those of the 

(non-projected) intrinsic states. 



In the top left of Fig. [T] is plotted the non-projected energy surface up to 7 
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MeV. The energy is normalized to the minimum we found and which turns out to 
be axial (/? = 0.26, 7 = 0.0°). In the top right of Fig. [T]we present the particle 
number restored (PNR) energy surface, where the projected energies are plotted at 
the deformation of the SCMF states they have been obtained fromQ One observes 
that the PNR does not change the global shape of the surface. In particular, the 
minimum remains axial. The PNR gives an almost constant overall energy gain of 
about 1 MeV for the entire surface. When projecting on both particle number and 
angular momentum, we get for J = |, the J which gives the absolute minimum 
after projection, the surface plotted in the bottom of Fig. [l] And as one can see, 
the topography of the surface has dramatically changed: the minimum has moved 
into triaxality (/3 = 0.33, 7 = 15.3°), and the surface is now centered around this 
minimum, which is about 5 MeV lower than the non-projected one. The angular- 
momentum projected surface is thus more rigid against oblate deformation than 
the non-projected one. 

In Fig. 2 we compare the J'^ = | energy surface of "^^Cr (left), now renormalized 
to its minimum, with the 2'" — 0^ energy surface of "^^Cr (right), also normalized to 
its minimum. One immediatly sees that the deepest part of the surface is very similar 
in both cases, in particular the minimum does not move much by the addition of 
one extra nucleon to ^®Cr. By contrast, the energy surface of the J'^ = | in "'^Cr 
is clearly more rigid than the J=0^ surface in "^^Cr, i.e. for the former the energy 
grows faster when moving away from the minimum. It will be interesting to see if 
these remarks remain true for different J coming from other blocked quasiparticles. 




0.0 0.2 0.4 0.6 0.0 0.2 0.4 0.6 

Fig. 2. Left: Energy surface, in the first sextant of tlie /3 — 7 plane, of tlie lowest J'^ = | after 
K mixing in *^Cr. Right: Same for the lowest J'^ = 0+ in **Cr. Energies are normalized to the 
minimum of each surface, which is marked by a filled circle. 



"^The same representation is also adopted for the angular-momentum restored surfaces. 
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Fig. 3. Top: square of the weights of the different J components summed over K, for three different 
deformations. Bottom: same but for the weights of the K components summed over J. 



4.2. J and K decompositions of SCMF states 

In order to better understand how the projection of quasiparticle states works, it 
is iUuminating to look at the decomposition of SCMF states on different angular 
momenta J and their intrinsic components K. We plot in Fig. [3] the square of the 
weights of J and K components for several points in the /3-7 plane with different 
deformation. If one looks first at the axially constrained point that corresponds to 
the minimum before projection (dark grey bars in Fig. [s]), one sees that there is 
indeed not just one J, but a full distribution beginning at J = |, which is the largest 
component in the distribution, and decreasing with increasing J. Together with the 
gain in energy, the fact that there are different J components in the SCMF state 
proves a posteriori the need for angular-momentum projection. One also notes that 
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there is no component with J = ^ or with J = |. This can be understood when 
looking at the K decomposition of this state. Here, one can see that there are only 
two different K components, K = | and its partneij^K = — |. As K has to be less 
than or equal to J, the J = 5 and J = | components vanish. It is also interesting 
to note that even if we don't have axial symmetry, a state constrained to axial 
deformation of its local density has in first approximation a good K. 
The light grey bars of Fig. [3] represent the J and K decomposition of the state 
which gives the triaxial minimum after projection on J = |. As one can see, the 
distribution has changed: it is now the J = | component which is the largest, and the 
distribution is more evenly divided around this maximum. But more importantly, 
the K decomposition one has now not only K = ± | , but a distribution of K with 
a staggering which favors K = ±i,±|,±|... So even if K = ± | remains the 
largest component in this distribution, this points out the necessity for a K mixing 
as presented in Eqns. Q and ([5]). 

Finally, the black bars in Fig. [3] represent the J and K decomposition of a state 
constrained to a spherical density distribution. One notices that there is only a 
J = I component in the decomposition. This can be easily understood: imposing 
spherical density, the single-particle wave functions are almost eigenstates of angular 
momentum, such that the HFB state is also almost an eigenstate of J, with the 
eigenvalue determined by the j of the blocked particle, which in this case is in the 
f-j/2 shell. By contrast, nothing in our calculation fixes which magnetic substate in 
the fi/2 shell is blocked. As a consequence, we find an arbitrary combination that 
depends mainly on the initialization of the blocked HFB calculation. 

5. Summary and Outlook 

We have reported first results obtained from a method that describes properties of 
odd-A nuclei by particle-number and angular-momentum projected configuration 
mixing based of self-consistently blocked one-quasiparticle states and using the 
Skyrme EDF. 

We focussed here on symmetry restoration in a regularized MR-EDF framework. 
A symmetry-restored GCM calculation along the same lines is underway and will 
be reported elsewhere.^ Having a method that allows to treat odd- A and even-even 
nuclei on the same footing provides a versatile tool for numerous studies of nuclear 
structure, such as the coupling of single-particle and shape degrees of freedom, the 
analysis of signatures for shell structure such as separation energies, g factors, spec- 
troscopic quadrupole moments, and their evolution with N and or the analysis 
of the interplay of pairing correlations and fluctuations in shape degrees of freedom 
for the odd-even mass staggering. 

''Finding components of equal weight for itK is a consequence of our choice of having eigenstates 
of the X signature operator; Rx = e~''^''^ . 
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